Scale-free avalanche dynamics in the stock market 
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Abstract 

Self-organized criticality has been claimed to play an important role in many natural and social systems. In the 
present work we empirically investigate the relevance of this theory to stock-market dynamics. Avalanches in 
stock-market indices are identified using a multi-scale wavelet-filtering analysis designed to remove Gaussian noise 
from the index. Here new methods are developed to identify the optimal filtering parameters which maximize the 
noise removal. The filtered time series is reconstructed and compared with the original time series. A statistical 
analysis of both high-frequency Nasdaq E-mini Futures and daily Dow Jones data is performed. The results of this 
new analysis confirm earlier results revealing a robust power law behaviour in the probability distribution function 
of the sizes, duration and laminar times between avalanches. This power law behavior holds the potential to be 
established as a stylized fact of stock market indices in general. While the memory process, implied by the power 
law distribution of the laminar times, is not consistent with classical models for self-organized criticality, we note 
that a power-law distribution of the laminar times cannot be used to rule out self-organized critical behaviour. 
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1. Introduction 

Attracted by several analogies with the dynam- 
ics of natural systems, physicists, especially dur- 
ing the last decade, have attempted to understand 
the mechanism behind stock-market dynamics by 
applying techniques and ideas developed in their 
respective fields [1] . 

In this context, possible connections between 
self-organized criticality (SOC) and the stock mar- 
ket, or economics in general, have been investigated 
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theoretically [2,3,4,5,6]. The theory of SOC, orig- 
inally proposed in the late eighty's by Bak, Tang 
and Wiesenfeld (BTW) [7] to explain the ubiquity 
of power laws in Nature, is claimed to be relevant 
in several different areas of physics as well as bio- 
logical and social sciences [4,8]. The key concept of 
SOC is that complex systems - i.e. systems consti- 
tuted by many non- linear interacting elements - al- 
though obeying different microscopic physics, may 
naturally evolve toward a critical state where, in 
analogy with physical systems near the phase tran- 
sition, they can be characterized by power laws. 
The critical state is an ensemble of metastable con- 
figurations and the system evolves from one to an- 
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other via an avalanche-like dynamics [8] . 

The classical example of a system exhibiting 
SOC behaviour is the 2D sandpile model [7,8,4]. 
Here the cells of a grid are randomly filled, by 
an external driver, with "sand" . When the gradi- 
ent between two adjacent cells exceeds a certain 
threshold a redistribution of the sand occurs, 
leading to more instabilities and further redistri- 
butions. The benchmark of this system, indeed of 
all systems exhibiting SOC, is that the distribu- 
tion of the avalanche sizes, their duration and the 
energy released, obey power laws. As such, they 
are scale-free. 

In the present work we search for imprints of 
SOC in the stock market by studying the statis- 
tics of the coherent periods (that is, periods of 
high volatility), or avalanches, which character- 
ize its evolution. We analyze the tick-by-tick be- 
haviour of the Nasdaq E-mini Futures (NQ) in- 
dex, P(t), from 21/6/1999 to 19/6/2002 for a to- 
tal of 2 19 data. In particular, we study the loga- 
rithmic returns of this index, which are defined as 
r(t) = ]n[P(t)/P(t-l)]. Possible differences be- 
tween daily and high frequency data have also been 
taken into consideration through the analysis of 
the Dow Jones daily closures (DJ) from 2/2/1939 
to 13/4/2004, for a total of 10 14 data. 

This work extends our earlier work on this sub- 
ject [9] by introducing new criteria to optimize the 
filtering of the time series essential to separating 
quiescent and avalanche dynamics. The properties 
of the time series reconstructed from the filtered 
returns are also examined. The issue regarding the 
presence of SOC in the stock market is of not only 
of theoretical importance, since it would lead to 
improvements in financial modeling, but could also 
enhance the predictive power [10] of Econophysics. 

In the next section we present the analysis 
methodology while in Sec. 3 the results of the anal- 
ysis are presented. Discussions and conclusions 
are contained in the last section. 



2. Avalanche Identification via Wavelets 

The logarithmic returns of stock indices rarely 
display intervals of genuinely quiescent periods, yet 



such periods are vital to the quantitative identifi- 
cation of avalanche dynamics. As such, noise must 
be filtered from the time series. Ideally, only Gaus- 
sian noise, associated with the efficient phases of 
the market where the movements can be well ap- 
proximated by a random walk [1] , is to be filtered 
from the time-series returns. Such dynamics have 
no memory and contrast the avalanche dynamics, 
i. e. anomalous periods characterized by large fluc- 
tuations, that we aim to analyze. 

Naively, one might simply set a threshold for 
the logarithmic returns, below which the index is 
deemed to be laminar. However, a simple threshold 
method is not appropriate, as it would include in 
the filtering some non-Gaussian returns at small 
scales that are relevant in our analysis. 

This difficulty is illustrated in Fig. 1 (Top) 
where the probability distribution function (PDF) 
for the returns of the NQ index, filtered using a 
fixed threshold of r t h — 5 standard deviations is 
shown by the open squares. In this case broad 
wings, related to events that do not follow Gaus- 
sian statistics, are clearly evident. 

However, an important stylized fact of financial 
returns - the intermittency of financial returns [1] 
- can be used to identify an appropriate filter- 
ing scheme. Already, physicists have drawn analo- 
gies with the well known phenomenon of intermit- 
tency in the spatial velocity fluctuations of hy- 
drodynamic flows [11,12,13]. Both systems display 
broad tails in the probability distribution func- 
tion [1] and a non-linear multifractal spectrum [12] 
as a result of this feature. The empirical analogies 
between turbulence and the stock market suggest 
the existence of a temporal information cascade for 
the latter [12] . This is equivalent to say that various 
traders require different information according to 
their specific strategies. In this way, different time 
scales become involved in the trading process. 

In the present work we use a wavelet method in 
order to study multi-scale market dynamics. The 
wavelet transform is a relatively new tool for the 
study of intermittent and multifractal signals [14]. 
This approach enables one to decompose the signal 
in terms of scale and time units and so to separate 
its coherent parts - i. e. the bursty periods related 
to the tails of the PDF - from the noise-like back- 
ground. This enables an independent study of the 
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avalanches and the quiescent intervals [15]. 

The wavelet transform (WT) is defined as the 
scalar product of the analyzed signal, f(t), at 
scale A and time t, with a real or complex "mother 
wavelet", tp(t). In the discrete wavelet transform 
(DWT) case, used herein, this reads: 

W T f{t) = ^=f f{u)i,(^^j du, (1) 

= 2-^ 2 J f(u) i)(2- j u -n)du, 

where the mother wavelet is scaled using a dyadic 
set. One chooses A = 2-?, for j = 0, L — 1, where 
A is the scale of the wavelet and L is the number of 
scales involved, and the temporal coefficients are 
separated by multiples of A for each dyadic scale, 
t = n2-?, with n being the index of the coefficient 
at the jth scale. 

The wavelet coefficients are a measure of the 
correlation between the original signal, f(t), and 
the mother wavelet, tp(t), at scale j and time n. 
In the analysis presented in the next section, we 
use the Daubechies-4 wavelet as the orthonormal 
basis [16]. However, tests performed with different 
sets do not show any qualitative difference in the 
results. 

The utility of the wavelet transform in the study 
of turbulent signals lies in the fact that the large 
amplitude wavelet coefficients are related to the ex- 
treme events corresponding to the tails of the PDF, 
while the laminar or quiescent periods are related 
to the coefficients with smaller amplitude [17]. In 
this way, it is possible to define a criterion whereby 
one can filter the time series of the coefficients de- 
pending on the specific needs. In our case, we adopt 
the method used in Ref. [17] and originally pro- 
posed by Katul et al. [18]. In this method wavelet 
coefficients that exceed a fixed threshold are set to 
zero, according to 

i otherwise . 

Here (. . .)„ denotes the average over the time pa- 
rameter n at a certain scale j and C is the threshold 
coefficient. In this way only the dynamics associ- 
ated with the efficient phases of the market where 
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Fig. 1. (Top) Comparison between the PDF of the orig- 
inal time series for the NQ index and its wavelet filtered 
version for C = 1. A Gaussian distribution is plotted for 
visual comparison. The simple threshold, r t ^ = 5, method 
for filtering is also shown. In this case it is clear that we do 
not remove just Gaussian noise, but also coherent events 
that can be relevant for the analysis. (Bottom) A window 
of the time series of the residuals obtained by subtracting 
filtered time series from the original time series. Avalanches 
of high volatility contrast periods of genuinely quiescent 
behavior. All the data in the plots have been standardized, 
r(t) — > (r(t) — (r))/a(r), where (...) and <r(r) are, respec- 
tively, the average and the standard deviation during the 
period under study. 

the movements can be well approximated by a ran- 
dom walk [1] are preserved. 

Once we have filtered the wavelet coefficients 
Wj iTl an inverse wavelet transform is performed, 
obtaining what should approximate Gaussian 
noise. The PDF of this filtered time scries is shown, 
along with the original PDF in Fig. 1 (Top). It is 
evident how the distribution of the filtered signal 
matches perfectly a Gaussian distribution. 

In the same figure (Bottom), we also show the 
logarithmic returns, R(t) , of the original time series 
after the filtered time series has been subtracted. 
Truly quiescent periods are now evident, contrast- 
ing the the bursty periods, or avalanches, which we 
aim to study. 

The time series of logarithmic prices is recon- 
structed from the residuals in Fig. 2 and is con- 
trasted with the one reconstructed from the filtered 
Gaussianly distributed returns. Note how, in the 
latter case, the time series is completely indepen- 
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Fig. 2. The Dow Jones time series is superimposed with 
the time series reconstructed from the filtered returns and 
the residual returns remaining after the filtered returns are 
subtracted from the original returns. The price behaviour 
generated by the "efficient" or filtered returns is largely 
independent of the observed price. The filtering parameter, 
in this case, is C = 1. 



Fig. 3. The excess of kurtosis plotted as a function of the 
filter parameter C for the NQ and DJ indices. A sample 
of Gaussian noise is also included for contrast. An optimal 
value of C ~ 1 is found, optimally filtering the original 
market time series to the level of noise. 

3. Data Analysis 



dent of the actual market price. 

To this point, the filtering parameter, C, has 
been constrained to 1, thus preserving coefficients 
that are less than the average coefficient at a par- 
ticular scale. However, one might wonder if it is 
possible to tune this parameter to maximally re- 
move the uninteresting Gaussian noise from the 
original signal. 

Fig. 3 illustrates the extent to which the filtered 
signal is Gaussian as a function of the filtering pa- 
rameter C. Here we report the value of the excess of 
kurtosis, K e = (r 4 ) /{r 2 ) 2 — 3, where (...) is the av- 
erage of the filtered time scries over the period con- 
sidered. For pure Gaussian noise this value should 
be 0. With this test we are able to identify C ~ 
1 as optimal for both the NQ and DJ indices in- 
vestigated here. Pure noise signals are completely 
filtered with this simple consideration: an exami- 
nation of the standard autocorrelation function of 
the filtered time series shows a complete tempo- 
ral independence, further confirming that we have 
successfully filtered Gaussian noise. 



Once we have isolated and removed noise from 
the time series we are able to perform a reliable 
statistical analysis on the avalanches of the residual 
returns, Fig. f (Bottom). In particular, we define 
an avalanche as the periods of the residual returns 
in which the volatility, v(t) = \r(t)\, is above a 
small threshold, typically two orders of magnitude 
smaller than the characteristic return. 

A parallel between avalanches in the classical 
sandpile models (BTW models) exhibiting SOC [7] 
and the previously defined coherent events in the 
stock market is straightforward. In order to test 
the relation between the two, we make use of some 
properties of the BTW models. In particular, we 
use the fact that the avalanche size distribution and 
the avalanche duration are distributed according 
to power laws, while the laminar, or waiting times 
between avalanches are exponentially distributed, 
reflecting the lack of any temporal correlation be- 
tween them [19,20]. This is equivalent to stating 
that the triggering process has no memory. 

Similar to the dissipated energy in a turbulent 
flow, we define an avalanche size, V, in the mar- 
ket context as the integrated value of the squared 
volatility over each coherent event of the residual 
returns. The duration, D, is defined as the inter- 
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Fig. 4. PDFs for the avalanche size, V, for the opti- 
mally-filtered NQ (Top) and DJ (Bottom) indices. The ex- 
ponents of the power laws are 7 ~ —2.4 (NQ) and 7 ~ — 1.9 
(DJ). 

val of time between the beginning and the end of 
a coherent event, while the laminar time, L, is the 
time elapsing between the end of an event and the 
beginning of the next. 

The results for the statistical analysis of the 
optimally-filtered NQ and DJ indices are shown in 
Figs. 4, 5 and 6 for the avalanche size, duration and 
laminar times, respectively. A power law relation is 
clearly evident for all three quantities investigated. 

The data analyzed herein display a distribution 
of laminar times different from the BTW model 
of the classical sandpile. As explained previously, 
the BTW model shows an exponential distribu- 
tion for L, derived from a Poisson process with no 
memory [19,20]. The power law distribution found 
here implies the existence of temporal correlations 
between coherent events. However this correlation 
may have its origin in the driver of the market, con- 
trasting the random driver of the classical sandpile. 



4. Discussion and Conclusion 

We have investigated the possible relations be- 
tween the theory of self-organized criticality and 
the stock market. The existence of a SOC state 
for the latter would be of great theoretical impor- 
tance, as this would impose constraints on the dy- 
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Fig. 5. PDFs for the avalanche duration, D, for the opti- 
mally-filtered NQ (Top) and DJ (Bottom) indices. The ex- 
ponents of the power laws are 7 ~ —4.2 (NQ) and 7 ~ —3.5 
(DJ). 




Fig. 6. PDFs for the laminar times between avalanches, 
L, for the optimally-filtered NQ (Top) and DJ (Bottom) 
indices. The exponents of the power laws are 7 ~ —2.1 
(NQ) and 7 ~ -2.3 (DJ). 

namics, as implied by the presence of a bounded 
attractor in the state space. Moreover, it would be 
possible to build new predictive schemes based on 
this framework [10]. 

After a multiscale wavelet filtering, an avalanche- 
like dynamics has been revealed in two samples 



of market data. The avalanches are characterized 
by a scale-free behaviour in the size, duration and 
laminar times. The power laws in the avalanche 
size and duration are a characteristic feature of a 
critical underlying dynamics in the system. 

However, the power law behavior in the lami- 
nar time distribution implies a memory process 
in the triggering driver that is absent in the clas- 
sical BTW models, where an exponential behav- 
ior is expected. Remarkably, the same features 
have been also observed in other physical con- 
texts [19,17,21,22,23]. 

The problem of temporal correlation in the 
avalanches of real systems, has raised debates in 
the physics community, questioning the practical 
applicability of the SOC framework [24]. Motived 
by this issue, several numerical studies have been 
devoted to including temporal correlations in SOC 
models [25,26,27]. A power-law distribution in the 
laminar times has been achieved, for example, 
by substituting the random driver with a chaotic 
one [28,29]. Alternatively, it has been shown that 
non-conservative systems, as for the case of the 
stock market, could be in a near-SOC state where 
dissipation induces temporal correlations in the 
avalanches while the power law dynamics persist 
for the size and duration [30,31]. 

In conclusion, a definitive relation between SOC 
theory and the stock market has not been found. 
Rather, we have shown that a memory process is 
related with periods of high activity. The memory 
could result from some kind of dissipation of infor- 
mation, similar to turbulence, or have its origin in a 
chaotic driver applied to the self-organized critical 
system. While a combination of the two processes 
can also be possible, it is the latter property that 
prevents one from ruling out the possibility that 
the stock market is indeed in a SOC state [29]. 

Similar power-law behaviour has been found in 
the ASX index for the Australian market [9] and 
different single stock time series. If this power- 
law behaviour is confirmed by further studies, this 
should be considered as a stylized fact of stock mar- 
ket dynamics. 
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